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Abstract 

A theory of exponential modified gravity which explains both early-time inflation and 
late-time acceleration, in a unified way, is proposed. The theory successfully passes the local 
tests and fulfills the cosmological bounds and, remarkably, the corresponding inflationary 
era is proven to be unstable. Numerical investigation of its late-time evolution leads to the 
conclusion that the corresponding dark energy epoch is not distinguishable from the one for 
the ACDM model. Several versions of this exponential gravity, sharing similar properties, 
are formulated. It is also shown that this theory is non-singular, being protected against the 
formation of finite-time future singularities. As a result, the corresponding future universe 
evolution asymptotically tends, in a smooth way, to de Sitter space, which turns out to be 
the final attractor of the system. 



I Introduction 

Modified gravity is getting a lot of attention from the scientific community, owing in particular to 
the remarkable fact that it is able to describe early-time inflation as well as the late-time (dark 
energy) acceleration epoch in a unified way. This approach appears to be very economical, as 
it avoids the introduction of any extra dark component (inflaton or dark energy of any kind) 
for the explanation of both inflationary epochs. Moreover, it may be expected that, with some 
additional effort, it will be able to provide a reasonable resolution of the dark matter problem 
as well as of reheating, two other important issues in the description of the evolution of our 
universe. Furthermore, as a by-product, modified gravity has the potential to lead to a number 
of interesting applications in high-energy physics. In particular, gravity, which provides a 
simple example of modified gravity, may serve for the unification of all fundamental interactions, 
including quantum gravity, in an asymptotically-free theory [IJ. Modified gravity may also be 
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used to provide the scenario for the resolution of the hierarchy problem of high energy physics [5]. 
Finally, the corresponding string/M-theory approach modifies gravity already in the low-energy 
effective-action approximation, so that a theory of the kind considered appears to be quite natural 
from very fundamental considerations (see, for instance, [3J). 

Presently, a number of viable F(R) gravities leading to a unified description as explained have 
been identified (for a recent review, see [4j, and for a description of the observable consequences 
of such models, see [5 ). It goes without saying that all those models are constrained to obey 
the known local tests, as well as cosmological bounds. However, this might not be such a severe 
problem, since already the first model proposed [6J which unified inflation with dark energy already 
satisfied many of these local tests. The real internal problem of F{R) gravity is related with its 
being a higher-derivative theory, which renders it highly non-trivial. This means that it is hard, 
in fact, to explicitly work with such theories and to get observable predictions from them. 

The main aim of this paper is to propose a reasonably simple but indeed viable version of F{R) 
gravity which consistently describes the unification of the inflationary epoch with the dark energy 
stage, while satisfying the known local tests and cosmological bounds. Speciflcally, to address the 
issues above, we here propose exponential gravity, which on top of being simple is moreover free 
from any kind of flnite-time future singularity and exhibits other very interesting properties, as 
we will see. 

The paper is organized as follows. In the next section we briefly review F{R) gravity as well 
as the corresponding FRW cosmological equations. Special attention is paid to de Sitter and 
spherically-symmetric solutions. Sect. Illll is devoted to the discussion of the viability conditions in 
F{R) gravity. These conditions are investigated for the simple and realistic theory of exponential 
gravity, proposed as a dark energy model, in Sect. IIVI In Sect.|V|we carry out a detailed analysis 
of our explicit proposal: exponential gravity which describes in a natural, unifying way both early- 
time inflation and late-time acceleration. It is there demonstrated, too, that the model leads to 
a satisfactory, graceful exit from inflation (the de Sitter inflationary solution being unstable). In 
Sect. IVII we show that the theory does not lead to any sort of finite-time future singularities. A 
careful numerical investigation of late-time cosmological dynamics is carried out in Sect. IVIII It 
will be demonstrated there that exponential gravity makes specific predictions which are not 
distinguishable from those of the ACDM model in the dark energy regime. The asymptotic 
behavior of the theory at late times is investigated in Sect. IVIIII Section IIXI is devoted to a 
somehow different model, a variant of exponential gravity which unifies unstable inflation with 
the dark energy epoch and which is protected against future singularities by construction. This 
opens the window to other variations of the basic model sharing all its good properties. In the 
discussion section |Xl a flnal summary and outlook are provided, and there is an Appendix on the 
Einstein frame. 



II gravity: general overview and FRW cosmology 



II. 1 The classical action 

The action of modified F{R) theories is j4]: 



5* = 



FjR) 
2k2 



^(matter) 



(III) 



where g is the determinant of the metric tensor g^^, £(mattcr) matter Lagrangian and F{R) 

a generic function of the Ricci scalar, R. In this paper we will use units where fee = c = /i = 1 
and denote the gravitational constant — SttGn = 87r/Mpj, with the Planck mass being Mp\ = 



G 



-1/2 



N 



1.2 X lO^^GeV. We shall write 



F{R) = R + f{R) . 



(II.2) 



The modification is represented by the function f{R) added to the classical term R of the Einstein- 
Hilbert action of General Relativity (GR) . In what follows we will analyze modified gravity in this 
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form, explicitly separating the contribution of GR from its modification. 

By variation of Eq. (jll.ip with respect to g^,y, we obtain the field equations: 



R,. - \r9,u - {rfj' + f . (II.3) 



Here, R^i, is the Ricci tensor and the part of modified gravity is formally included into the 'modified 



gravity' stress-energy tensor T^*^, given by 



TfJ" = ^^J^ \]^9^^AF{R) - RF\R)] + (V^V. - g^,U)F'{R)^ 



(II.4) 



The prime denotes derivative with respect to the curvature i?, is the covariant derivative 
operator associated with g^^ and Dc^ = g'^^'V^V^^ is the d'Alembertian for a scalar field (j). 
^^(mattci) gj^^gjj ]-,y |-]^g non- minimal coupling of the ordinary matter stress-energy tensor 
with geometry, namely, 

/T^(mattcr) /Ti(mattcr) /tt r\ 



In general, j'^™^"''') = diag(p,p,p,p), where p and p are, respectively, the matter energy-density 
and pressure. When F(i?) = i?, T^G ^ q and f^'°-"") = y^-attor)^ 

It should be noted that, due to the diffeomorphism invariance of the total action, only 
is covariantly conserved and, formally, -pij]^ may be interpreted as an effective gravitational 
constant, assuming we are dealing with models such that F'{R) > 0. 

The trace of Eq. (|II.3|) reads 

3aF'{R) + RF'{R) - 2F{R) = K2j,(matter) ^ ^jj 

with 7"(inattor) |-]^g trace of the matter stress-energy tensor. We can rewrite this equation as 



where 

dV,s ^ 1 
dF'{R) 3 



2F(i?) - + K^T^™^"")] , (II 



F'{R) being the so-called 'scalaron' or the effective scalar degree of freedom. On the critical points 
of the theory, the effective potential Vcs has a maximum (or minimum), so that 

□F'(i?cp) = 0, (II.9) 

and 

2F(i?cp) - RcpF'iRcp) = -K^T^'"""'^'-) . (11.10) 

Here, Rcp is the curvature of the critical point. For example, in absence of matter, i.e. y(niattcr) _ 
0, one has the de Sitter critical point associated with a constant scalar curvature i?dSj such that 

2F(i?ds) - RdsF'iRds) - . (n.ll) 

Performing the variation of Eq. pi.6l) with respect to R, by evaluating \I\F'{R) as 

DF'{R) = F"{R)DR + F"'V^'RV^R , (11.12) 
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we find, to first order in SR, 



1 



F"{R) ■ 



-nsR + 1 

F"'{R) 



F"'{R) F"'{R) 



F"{R) \ F"{R) 



3F"{R) 

2" 



[2F{R) - RF'{R) + K^T™^""] 
R F'{R) 



g'^'^V^RV.R 



+2 



3(F"(i?))2 
F"'iR) 



[2F{R) - RF'{R) + K^T™^""'] 



3 3F"{R) 

^2 ^ymatter 
3F"{R) dR 



SR 



F"{R) 



g^'^V^.RVJR + 0{SR^) ~ . 



(11.13) 



This equation can be used to study perturbations around critical points. By assuming R — Rq 
const (local approximation), and SR/Rq <^ 1, we get 



aSR ~ m^SR + 0{5R^) 



where 



F"iRo) 



F"{Ro) dR 



Ro 



Here, Eq. pi.lOp with Rqp — Rq has been used. Note that 

2 SVeff 



dF'{Rf 



(11.14) 
(11.15) 

(11.16) 



The second derivative of the effective potential represents the effective mass of the scalaron. Thus, 
if m? > (in the sense of the quantum theory, the scalaron, which is a new scalar degree of 
freedom, is not a tachyon), one gets a stable solution. For the case of the de Sitter solution, 
is positive provided 

^'^^"^^ > 1 (II 17) 

The value of the above mass will be later used to check the emergence of the Newtonian regime. 

II. 2 Modified FRW dynamics 

The spatially-flat Friedman-Robertson- Walker (FRW) space-time is described by the metric 

ds^ = -dt^ + a'^{t)d^^ , (11.18) 
where a{t) is the scale factor of the universe. The Ricci scalar is 

i? = 6 (2H^ + H'j . (11.19) 

In the FRW background, from (/i, ly) — (0, 0) and the trace part of the (/i, i^) = {i,j) = 1, 3) 



(11.20) 



components in Eq. pi.3p . we obtain the equations of motion: 



Poff 



4 ( 2i? + 



(11.21) 



where pcff and Pcff are the total effective energy density and pressure of matter and geometry, 
respectively, given by 



Pcff 



F'iR) 



2k2 



{F'{R)R - F{R)) - 6HF'{R) 



(11.22) 



4 



Pcff 



p 



1 

2^ 



-{F'{R)R - F{R)) + AHF'{R) + 2F'{R) 



Here, H = a{t)/a(t) is the Hubble parameter and the dot denotes thue derivative d/dt. 
the form that the total stress tensor in Eq. (jH.3[) assumes in the FRW space-time. 
The standard matter conservation law is 



(n.23) 
This is 



and, for a perfect fluid. 



uj being the thermodynamical EoS-parameter of matter. 

We also introduce the effective EoS by using the corresponding parameter Woff 



Wcff 



PeS 
PeS 



and get 



Woff 



-f - 



2H 
3IP 



If the strong energy condition (SEC) is satisfied (woff > —1/3), the universe expands in a 
ated way, and vice-versa. We are interested in the accelerating FRW cosmology below. 



(11.24) 
(n.25) 

(n.26) 

(n.27) 

deceler- 



II. 3 Spherically symmetric solutions 

In this section we investigate spherically-symmetric solutions (like the Schwarschild black hole), 
which constitute an essential element for the local tests of modified gravity under consideration. 
For the metric, we start from a static, spherically symmetric ansatz of the type. 



-C(r)e2"('')dt2 



CV) 



r'^dn . 



where dil — r^{d9^ + sin^ OdcjP) and a{r) and C(r) are functions of r. 
Plugging this ansatz into the action (|II.ip . and noting that 



(n.28) 



ar 



—a{r) - 2C(r) 
dr 



dr 



- C{r)-2C{r)—a{r) 
dr'' dr'^ 



^ d r( ^ 4C(r) d 2C{r) 

— -rC [r) —a [r) — 

r dr r dr 



one arrives at the following equations of motion (in vacuum) [7]: 



RspF'{Rsp) — F(i?sp) 
F'(Rsp) 
2C(r)f "(i?,p) IfR^ 
F'(i?sp) dr^ 



(1 - C(r) - r{dC{r)/dr)) 



2 dC{r)/dr \ dR^ 
r ^ 2C(r) I 



F"'(i?,p) fdR, 



F"{Rsp) V dr 



(11.29) 



(n.30) 



0, 



dajr) /2 F"{R,p)dR,p 
dr \r F'{Rsp) dr 



F"(i?sp) d^Rsp F'"{R,p) fdR 



sp 



F'iRsp) dr^ F'iRsp) V dr 



0. 



(H.ai) 



These equations form a system of ordinary differential equations in the three unknown quantities 
a(r), C{r), and Rsp{r). When F{R) = R, the above system of differential equations lead to the 
Schwarzschild solution, namely 

a{r) — const , (11.32) 
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C{r) =(^l-^y (11.33) 

with M a dimensional constant, and Rgp = 0. 

Another well known case is the one associated with i?sp being constant. As a result, with 
a — 0, Eq. (jll.Sip is trivially satisfied, and the other two equations lead to the Schwarschild-de 
Sitter solution 

cM.(i-HM_^), ,„,4, 

with 

-Rsp = 4A , (11.35) 

and 

2F(i?sp) = RspF'iRsp) . (11.36) 



III Viability conditions in F(i?)-gravity 

The viability conditions follow from the fact that the theory is consistent with the results of 
General Relativity if F(0) = 0. In this way we can have the Minkowski solution. Recall that 
in order to avoid anti-gravity effects, it is required that F'{R) > 0, namely, the positivity of the 
effective gravitational coupling. 



III.l Existence of a matter era and stability of cosmological perturba- 
tions 

On the critical points, F'{R) = 0, and from Eqs. pi.22p and pi.23|) . one has 

^ F^) + 2^ [(^'(^)^ - ^(^))]} ' 

P^*^ " + i [-(^'(^)^ - Pirn] ■ (III-2) 

During the matter dominance era, we havepeff = 0. As as consequence, neglecting the contribution 
of the radiation, namely p = 0, one has peS = p/ F'{R) and 

(III.3) 



thus 



and using Eq. (|nL3|) . this leads to 



F{R) 

d_ r RF'jR) 
dR V F{R) 

F"{R) 



0, (III.4) 



= 0, (III.5) 



F'(i?) 

so that, during the matter era, we have F"{R) ~ 0. 

In order to reproduce the results of the standard model, where R = p when matter drives the 
cosmological expansion, a F(i?)-theory is acceptable if the modified gravity contribution vanishes 
during this era and F'{R) ~ 1. However, another condition is required on the second derivative 
of F{R): it not only has to be very small, but also positive. This last condition arises from the 
stability of the cosmological perturbations. We consider a small region of space-time in the weak- 
field regime, so that the curvature is approximated by i? ~ i?o + SR. where i?o — const. From 
Eq. (jlLlOp . we obtain 

- K^T = F'{R)R + 2{F{R) - RF'{R)) (III.6) 
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and, since F'{R) ~ 1 and \ {F{R) — R)\ ^ R, we can expand this expression as 

- K^T ~= -K^ {T\r + dT) = R + 2{F'{R) - l)SR , (III.7) 
with (SR/R) <C 1. By evaluating it at i? = i?o, from Eq. (III.15|) one has 

2^,1 f 2 F'iRo) ^\ ^,1 f 1 (F'iRo) - 1) 



^ 3 [f"{Ro) F"(i?o) ^ 3 [f"{Ro) F"(i?o) ) ' ^"^'^^ 

which is in agreement with Ref. [8J. Since {F'{Rq) — 1) <C 1 and as F"{Rq) is very close to zero, 
if F"{Rq) < the theory will be strongly mistable. Thus, we have to require F"{R) > during 
the matter era. 

111. 2 Existence and stability of a late-time de Sitter point 

It is convenient to introduce the following function, G{R), 

G{R) = 2F{R) - RF'{R) . (III.9) 

On the zeros of G{R) we recover the condition in Eq. (jll.lll) and we have the de Sitter solution 
which describes the accelerated expansion of our universe. If the condition in Eq. pi.l7[) is satisfied, 
the solution will be stable. 

A reasonable theory of modified gravity which reproduces the current acceleration of the uni- 
verse needs to lead to an accelerating solution for R = 4A, A being the cosmological constant 
(typically A ~ lO^^^eV^). Recall that in the de Sitter case, the EoS parameter ujcs — i-^ds = 
and all available cosmological data confirm that its value is actually very close to —1. The pos- 
sibility of an effective quintessence/phantom dark energy is not excluded, but the most realistic 
solution for our current universe is a (asymptotically) stable de Sitter solution. 

111. 3 Local tests and the stability on a planet's surface 

GR was first confirmed by accurate local tests at the level of the Solar System. A theory of 
modified gravity has to admit an asymptotically flat (this is important in order to define the mass 
term) static spherically-symmetric solution of the type pi. 341) . with A very small. The typical 
value of the curvature in the Solar System, far from sources, is R = R*, where R* ~ 10~^"'^eV^ 
(it corresponds to one hydrogen atom per cubic centimeter). If a Schwarzshild-de Sitter solution 
exists, it will be stable provided 

-I^i^>l. (III.IO) 

R*F"{R*) ^ ' 

The stability of the solution is necessary in order to find the post- Newtonian parameters in GR. 

Concerning the matter instability [HI [SI IS] j this might also occur when the curvature is rather 
large, as on a planet (i? ~ 10~'^^eV^), as compared with the average curvature of the universe 
today (i? ~ 10~^^eV^). In order to arrive to a stability condition, we can start from Eq. (III.13|) . 
where R = Rb assumes the typical curvature value on the planet and SR is a perturbation due to 
the curvature difference between the internal and the external solution. Since 
and SR depends on time only, one has 



dt{SR) ~ const + U{Rb)SR , (III.ll) 



where 




F"'{Rb) 
F"{Rb) 



g'^^S/f.RbV^Rb - 



Rb , F'(Rb 



3 3F"{Rb) 



3{F"{RbW ' ^'^'^^'^ ' ^'^} ■ ^"^-^^^ 
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If U{Rb) is negative, then the perturbation SR becomes exponentially large and the whole system 
becomes unstable. Thus, the matter stability condition is 

U{Rb) > 0, where Rb ~ IQ-^^eV^ . (III.13) 

At the cosmological level this means that F"{R) ~ 0+, in the matter era. If F{R ~ Rb) ~ R, 
Eq. (|TlLT2ll reads, simply, U{Rb) = l/(3F"(i?h)). 



III. 4 Existence of an early-time acceleration and the future singularity 
problem 

In order to reproduce the early-time acceleration of our universe, namely the inflation epoch, the 
modified gravity models have to admit a solution for Woff in Eq. pi.27p smaller than —1/3. An 
important point is that this solution should be unstable. 

If the model reproduces the de Sitter solution when i? = i?ds - lO^o-^^GeV^ (th is is the 
typical curvature value at inflation), we have to require that Eq. (jII.17p is violated. Thus, the 
characteristic time of the instability ti is given by the inverse of the mass of the scalaron in 
Eq. (|lLT5l) : 



1 




m 





F'iRds) 



F'{Rds)-RdsF"iR. 



dS) 



(III.14) 



Note that if the scalaron mass is equal to zero, a more detailed analysis, as in Sect. IIX[ is needed. 

Furthermore, it is well- know that many of the effective quintessence/phantom dark energy 
models, including modified gravity, bring the future universe evolution to a finite-time singularity. 
The most familiar of them is the famous Big Rip [TU] , which is caused by phantom dark energy. 
Finite-time future singularities in modified gravity have been studied in Refs. |12[ [T5] . As is 
known, a finite-time future singularity occurs when some physical quantity (as, for instance, the 
scale factor, effective energy- density or pressure of the universe or, more simply, some of the 
components of the Riemann tensor) diverges. The classification of the (four) finite-time future 
singularities has been done in Ref. [TT] . Some of these future singularities are softer than others 
and not all physical quantities necessarily diverge on the singularity. 

The presence of a finite-time future singularity may cause serious problems to the cosmological 
evolution or to the corresponding black hole or stellar astrophysics [14 1. Thus, it is always necessary 
to avoid such scenario in realistic models of modified gravity. It is remarkable that modified 
gravity actually provides a very natural way to cure such singularities by adding, for instance, an 
i?^-term [151 112| . Simultaneously with the removal of any possible future singularity, the addition 
of this term supports the early-time infiation caused by modified gravity. Remarkably, even in 
the case inflation were not an element of the alternative gravity dark energy model considered, 
it will eventually occur after adding such higher-derivative term. Hence, the removal of future 
singularities is a natural prescription for the unified description of the inflationary and dark energy 
epochs. 



IV Realistic exponential gravity 

In Refs. [ini [TTl [TSj several versions of viable modified gravity have been proposed, the so-called 
one-step models, which reproduce the current acceleration of the universe. They incorporate a 
vanishing (or fast decreasing) cosmological constant in the flat [R — > 0) limit, and exhibit a 
suitable, constant asymptotic behavior for large values of R. The simplest one was proposed in 
Ref. ^ 

F{R) = R-2A(^1- e"-^/^") . (IV.l) 

Here, A ~ lO^^^eV'^ is the cosmological constant and i?o — A a curvature parameter. In flat space 
one has F(Q) — and recovers the Minkowski solution. For i? ^ i?o, F{R) ~ R — 2A, and the 
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theory mimics the ACDM model. Note that late-time cosmology of such exponential gravity was 
also considered in Ref. [19]. 
For simplicity, we will set 

f{R) = -2A(1 - e"^/^" ) , (IV.2) 

and thus 

fiR) = -2 Ae-«/«c , (IV.3) 

f"{R) = 2 Ac-«/«o (IV.4) 
Rq 

Since \f'{R ^ Rq)\ ^ 1, the model is protected against anti-gravity during the whole cosmological 
evolution, until the de Sitter solution (i? = 4A) of today's universe is reached. For large values of 
the curvature, F{R ^ Rq) ~ R, and we can reconstruct the matter-dominated era, as in GR. In 
particular, F"{R > Rq) ~ Ae"-^/^" ~ 0+, and we do not have any instability problems related to 
the matter epoch, obtaining matter stability on a planet's surface and at the solar system scale. 
Let us consider the G{R) function of Eq. pil.9|) . 

G{R) = R + 2f{R)-Rf'{R). (IV.5) 

Since G(0) = 0, one has a trivial de Sitter solution for R = 0. Consider now 

G'{R) = l + f{R)-Rr{R). (IV.6) 

Since G"(0) < 0, the function G{R) becomes negative and starts to increase after R = Rq. For 
R = 4A, F{R) ~ -2A, F'{R) ~ 1 and F"{R) ~ 0+. It means that G(4A) ~ and we find the de 
Sitter solution of the dark energy phase which is able to describe the current acceleration of our 
universe. After this stage, G{R > 4A) ~ i? > and we do not find other de Sitter solutions. Note 
that the de Sitter solution for R = 4A is stable, since the first term in Eq. (|II.17|) diverges. On the 
other hand, the Minkowski space solution is unstable. Summing up, we have two FRW-vacuum 
solutions, which correspond to the trivial de Sitter point for R = and to the stable de Sitter 
point of current acceleration, for R = 4A. 

Finally, we have to consider the existence of spherically-symmetric solution. In i? = we find 
the Schwarzschild solution, which is unstable. On the other hand, the physical Schwarzschild-de 
Sitter solutions are obtained for R ^ Rq. For example, in the Solar System, R* ~ 10~^^eV^. In 
this case {F{R*) ~ i?* — 2A) we find the Schwarzschild-de Sitter solution as in Ea. (|lL34ll . which 
can be approximated with the Schwarzschild solution of Eqs. (|TL32])-(|TL33|, owing to the fact that 
A is very small. For R* ^ Rq, Eq. (jllLlOp with Rsp — R* is satisfied and the solution is stable. 

The description of the cosmological evolution in exponential gravity has been carefully studied 
in Refs. [191 where it was explicitly demonstrated that the late-time cosmic acceleration fol- 
lowing the matter-dominated stage, as final attractor of the universe, can indeed be realized. By 
carefully fitting the value of Rq, the correct value of the rate between matter and dark energy of 
the current universe follows (see Sect. IVI[) . Our next step will be to generalize the model in order 
to describe inflation. We will follow the method first suggested in Ref. |18) . 



V Inflation 

A simple modification of the one-step model which incorporates the inflationary era is given by 
a combination of the function discussed above with another one-step function reproducing the 
cosmological constant during inflation. A quite natural possibility is 



=i?-2A(^l-e"Wj -Ai l^l-e V«J j+-fR". (V.l) 

For simplicity, we call 

/, --A Jl-e'(*)"V (V.2) 
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where Ri and A.^ assume the typical values of the curvature and expected cosmological constant 
durmg mflation, namely Ri, Ai ~ 10^'^~'^*eV^, while n is a natural number larger than one. The 
presence of this additional parameter is motivated by the necessity to avoid the effects of inflation 
during the matter era, when i? <C so that, for n > 1, one gets 

on 

(V.3) 

The last term in Eq. (jV.ip . namely 7i?", where 7 is a positive dimensional constant and a a real 
number, is necessary to obtain the exit from inflation. If 7 ^ and a > 1, the effects of 

this term vanish in the small curvature regime, when i? <C -Ri and 

Note that /^(O) = and /^(i? > ~ -A^. We also obtain 

fUR) = - \n e , (V.5) 

The first derivative //(i?) has a minimum a.t R ~ R, where fi{R) — 0. One gets 

Thus, in order to avoid the anti-gravity effects < 1), it is sufficient to require |/j'(-R)| < 1- 

This leads to 

R, > A,n f " e-^ . (V.8) 



For example, one can choose n ~ 4. In this case Eq. (jV.Sp is satisfied for i?^ > 1.522 A^. A 



reasonable choice is Ri — 2Ai. The last power-term of Eq. (jV.l[) does not give any problem with 
anti-gravity, because its first derivative is positive. 

It is necessary that the modification of gravity describing inflation does not have any influence 
on the stability of the matter era in the small curvature range. When i? <C the second 
derivative of such modification, namely 



n'{R)+a{a~l)jR^-'c,^ 



-n(n-l)(- 



(V.9) 



must be positive, that is 

n>a. (V.IO) 

We require the existence of the de Sitter critical point i?ds which describes inflation in the 
high-curvature regime of fi{R), so that fi{RdS ^ Ri) — —-^i and fl{Rds ^ Ri) — 0+. In this 



region, the role of the first term of Eq. (jV.l|l is negligible, while the term jR°' needs to be taken 
into account. For simplicity, we shall assume that 7 — ^/R^s^- The function G{R) in Eq. pil.9|) . 

G(i?)-i? + 2/,-i?/; + ^|=;^i?«, (v.ii) 

has to be zero on the de Sitter solution. We get 

2A- 

Rds = ir-^ , Rds » R^ ■ (V.12) 

6 — a 
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Since Ai ^ Ri, in order to satisfy the last two conditions simultaneously, one has to choose 



2<a<3. (V.13) 
Let us consider the effective scalar mass of Eq. pi.lSp on the de Sitter solution: 

2 i?ds 1 + 2a - q;2 ~ 



a{a — 1) 



(V.14) 



It is negative if a > 2.414. In this case inflation is strongly unstable. Using Eq. (IIII.14|) we derive 
the characteristic time of the instability as 

U 7^ 10"^" - 10"^^ sec , (V.15) 

V^dS 

in accordance with the expected value. The new condition on a, in order to have unstable inflation, 
is 

5/2<a<3. (V.16) 
Now, we will try to reconstruct the evolution of the function G{R) in Eq. (jV.llI) . 

G{R) = R + 2fdR) - RfliR) + (2 - . (V.17) 

When i? = 0, we find a trivial de Sitter point and G'(O) = 0. For the first derivative of G{R), 

G\R) = 1 + n{R) - Rfl'iR) + a{2 - a)^ . (V.18) 

^dS 

G'(0) > and G{R) increases. Since f,l'{R) starts being positive for R > R (where R is expressed 



as in Eq. (jV.7P ') and 2 — a < 0, it is easy to see that G{R) begins to decrease at around R = Ri 
and that it is zero when R = i?dS- After this point, G'{R > Rds) < and we do not have other de 
Sitter solutions. On the other hand, it is possible to have a fluctuation of G{R) along the i?-axis 
just before the de Sitter point describing inflation takes over. In order to avoid other de Sitter 
solutions (i.e., possible final attractors for the system), we need to verify the fulfillment of the 
following condition: 

G{R) > for < i? < i?ds ■ (V.19) 

Precise analysis of this condition leads to a transcendental equation. In the next subsection we 
will limit ourselves to a graphical evaluation. In general, it will be enough to choose n sufficiently 
large in order to avoid such effects. 



V.l Construction of a realistic model for inflation 

By taking into account all the conditions met in the previous paragraph, the simplest choice of 
parameters to introduce in the function of Eq. (jV.l[) is 

n = 4, a = ^, (V.20) 

while the curvature Ri is set as 

Ri = 2A,; . (V.21) 

In this way, n > a and we avoid undesirable instability effects in the small-curvature regime. Ri 
satisfies Eq. (jV.8|) and we have no anti-gravity effects. From Eq. (jV.12[) one recovers the unstable 
de Sitter solution describing inflation as 

Rds = 4A, . (V.22) 
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We note that, due to the large value of rt, Bds is sufficiently large with respect to Ri, and fi{Rds) — 
—Ai. One can also expect that, on top of this graceful exit from inflation, the effective scalar degree 
of freedom may also give rise to reheating, in analogy with Ref. [21\. 

In Fig. [T]a plot of G{R) is shown. The zeros of G{R) correspond to de Sitter solutions. One 
can see that the only non-trivial zero is the de Sitter point of Eq. (|V.22[) . and here the function 
crosses the i?-axis up-down, according to the instability of such solution (since F"{R > R) > 0, we 
get G'{R) ^ m^<0). This means that the inflationary de Sitter point corresponds to a maximum 
of the theory (without matter/radiation). The system gives rise to the de Sitter solution where the 
universe expands in an accelerating way but, suddenly, it exits from inflation and tends towards 
the minimal attractor at i? = 0, unless the theory develops a singularity solutions for i? — )• oo. 
In such case, the model could exit from inflation and move in the wrong direction, where the 
curvature would grow up and diverge, and a singularity would appear. In the next section we 
will recall some important facts about singularities, which will be considered in the context of 
exponential gravity. 

G 




Figure 1: Plot of G{R/Ai) The zeros correspond to the de Sitter solutions. 



VI Finite-time future singularities 

In general, future singularities appear when the Hubble parameter has the form 

where h and Iq are positive constants, and t < to, because it should correspond to an expanding 
universe. Here /3 is a positive constant or a negative non-integer number, so that, when t is close 
to H or some derivatives of _ff, and therefore the curvature, become singular. 

Note that such a choice of the Hubble parameter corresponds to an accelerated universe, 
because on the singular solution of Eq. (jVI.ip it is easy to see that the strong energy condition 
(Pcff + 3pcff > 0) is always violated when /3 > 0, or for small values of t when f3 < 0. What 
means that, in any case, a singularity could emerge at the flnal evolution stage of an accelerating 
universe. 

The finite-time future singularities can be classified as follows |11| : 

• Type I (Big Rip): for t — !• to, a{t) — oo, pcs oo and |poff| — > oo. It corresponds to f3 — 1 
and (3 > 1. 

• Type II (sudden): for t — > to, a{t) — >■ oq, PcS — >■ Pn and |poff| — > oo. It corresponds to 
-1< /3 < 0. 

• Type HI: for t — > tg, a{t) oq, PoS — ^ oo and |poff | — > oo. It corresponds to < /3 < 1. 



12 



• Type IV: for t — > to, a(0 0,0, PcS ^ 0, |pcff| and higher derivatives of H diverge. The 
case in which p and/or p tend to finite values is also included. It corresponds to /3 < — 1 but 
/? is not any integer number. 

Here, 00(7^ 0) and po 8^re positive constants. 

It is easy to see that in the case of Type I and Type III singularities we have R — > +00. Those 
types are the most dangerous in the cosmological scenario. On the other hand, also the soft Type 
II and Type IV {R — const) singularities may cause various problems related , for instance, to 
the description of stellar astrophysics. 



VI. 1 Singularities in exponential gravity 

Let us now consider the exponential model of Eq. (jV.l[) . but avoiding the last term, 7_R". When 
R +00, we have F{R — > +00) ~ i? + const, F'{R — >■ +00) ~ 1, while the high-order derivatives 
of F{R) tend to zero in an exponential way. This means that PeS in Eq. pi.22p . by neglecting the 
matter contribution, tends to a constant on the singular solution of Eq. (jVI.ip with /3 > (namely, 
R —7- -|-(xi). For this reason, neither Type I nor Type III singularities, where p^g has to diverge, 
can appear in exponential gravity. However, a model of this kind, which mimics the cosmological 
constant in the high-curvature regime, can be affected by Type II or IV singularities. 

Regarding Type II singularities, we have to consider the behavior of the model when R is 
negative, which is very different from the case of positive values of the curvature. This is the reason 
why this kind of singularities does not appear: when R — > —00, p^s of Eq. pi.22p exponentially 
diverge, and the sudden singularity, where pcff has to tend to a constant, cannot be realized. 

Also Type IV singularities are not realized: when R 0", F{R 0~) ~ i? — R/Rq, and it is 
easy to see that poff of Eq. (III.22P behaves as l/{to—t)^^^, and it is larger than H^{— /i^/(to — 1)^^), 
when /3 < —1. For this reason, Eq. (|II.20p is inconsistent with Type IV singularities. The argument 
is valid also if we consider the more general case where H = Hq + h/ {to — t)^ and tends to the 
positive constant Hq in the asymptotic singular limit: also in this situation F{R) approaches a 
constant like l/(to — l)*^^^, while the time dependent part of behaves as l/(to ^ t)^ ■ 

Consider now adding back the term jR". This becomes relevant just when R ^ Ri, so that 
it can produce Type I or HI singularities, only. However, in Ref. |13j it is explicitly demonstrated 
that the model R + jR" is protected against singularities of this kind if 2 > a > 1. 

Thus, we have found our theory to be free from singularities. In particular, when our model is 
approximated as 

F{R) ~ i? - A, + 7i?" , (VI.2) 

Type I or HI singularities do not occur. When inflation ends, the model moves to the attractor 
de Sitter point. In this way, the small curvature regime arises, the first term of Eq. (|V.ip becomes 
dominant and the physics of the ACDM model are reproduced. 



VII Dark energy epoch 

We will now be interested in the cosmological evolution of the dark energy density pde = 
PcS — p/F'iR) in the case of the two-step model of Eq. (IV. II) . near the late-time acceleration 
era describing the current universe. We assume the spatially-flat FRW metric of Eq. (III.lSp . Let 
us follow the method first suggested in Ref. [TB] and more recently used in Ref. [20j . 
To this end, we introduce the variable 

_ Pde H'^ _3 _4 

Here, is the energy density of matter at present time, is the mass scale 

m2 = ~ 1.5 X lO-e^eV^ , (VII.2) 
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and X is defined as 

(0) 

X EE 3.1 X 10-4^ (VII.3) 

Pra 

where p'r^ is the energy density of radiation at present (the contribution from radiation is also 
taken into consideration). 

The EoS-paraineter wde for dark energy is 

3 yn d(lna) 

By combining Eq. (|II.20p with Eq. (III. 191) and using Eq. (jVILip . one gets 

■h-TT, r + -'^2yH + ^3 = 0, (VII. 5) 



where 



d(lna)2 'd(lna) 

T = 1 2-F'{R) 

yn + a-s+Xa-^ 3m^F" (R) ' ^ ' 

y^ + a-^ + Xa-^ &fh^F"{R) ' 



and thus, we have 



i? = 3mM ^ + 4yH + ) . (VII.9) 



The parameters of Eq. (|V.1I) are chosen as follows: 

A= (7.93)TO^ 

A, = 101ooA, 
i?, = 2A, , n = 4 , 
_ 5 _ 1 

"~ 2 ' ~ (4A,)"-i ' 

i?o = 0.6A, 0.8A, A. (VII.IO) 

Eq. (jVII.Sp can be solved in a numerical way, in the range oi Rq ^ R Ri (matter era/current 
acceleration), yn is then found as a function of the red shift z, 

z=--l. (VII.ll) 
a 

In solving Eq. (jVII.5p numerically we have taken the following initial conditions at z = 

dyn 



d{z) 



= 0, 
A 



yn (VII.12) 

which correspond to the ones of the ACDM model. This choice obeys to the fact that in the high 
red shift regime the exponential model is very close to the ACDM Model. The values of Zi have 
been chosen so that RF"{z — Zi) ^ 10~^, assuming R — 'ifh?{z + l)'^. We have z,; = 1.5, 2.2, 2.5 
for Rq = 0.6A, 0.8A, A, respectively. In setting the parameters, we have used the last results of 
the WMAP, BAG and SN surveys [53]. 
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Figure 2: Plot of wde for i?o 0.6A. 
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Figure 3: Plot of wde for i?o 


= 0.8A. 
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Figure 4: Plot of wde for i?o = A. 



Using Eq. (|VII.4p . one derives wde from ?/h- In Figs. [5J [3] and HI we plot wde as a function 
of the redshift z for Rq — 0.6A, 0.8A, A, respectively. Note that wde is very close to minus one. 
In the present universe {z = 0), one has wde = -0.994, -0.975, -0.950 for i?o = 0.6A, 0.8A, 
A. The smaller i?o is, our model becomes more indistinguishable from the ACDM model, where 
cjde = -1- 
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We can also extrapolate the behavior of the density parameter of dark energy, SIde, 



^^DE — 



Pde _ m 

PcB yn + iz + lf + x{z + lY 



(VII.13) 



Plots of r^DE as a function of the redshift z for Bq = 0.6A, 0.8A, A, are shown in Figs. [5j Eland 
in For the present universe (z = 0), one has Que = 0.726, 0.728, 0.732 for Rq = 0.6A, 0.8A, A, 
respectively. 
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Figure 5: Plot of floE for Rq = 0.6A. 
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Figure 6: 


Plot of Ode for Rq = 0.8A. 



The data are in accordance with the last and very accurate observations of our present universe, 
where: 



I^DE - -0.972^0.060 ' 

f^DE = 0.721 ±0.015. 



(VII.14) 



As last point, we want to analyze the behavior of the Ricci scalar in Eq. (|VII.9[) for Rq — 0.6A, 
0.8A, A. Results are shown in Figs. \E\ and [TUl We clearly see that the transition crossing the 
phantom divide does not cause any serious problem to the accuracy of the cosmological evolution 
arising from our model. In particular, R{z — > —1^) tends to 12rh^yii{z —1^), which is an 
effective cosmological constant (note that Rq is small and we are close to the value of the ACDM 
model, where 12rh?y^ = 4A). As a consequence, the de Sitter solution is a final attractor of our 
system and describes an eternal accelerating expansion. 
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Figure 7: Plot of JIde for Rq = A. 
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Figure 8: Plot of R/A for Rq = 0.6A. 
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Figure 9: Plot of R/A for Rq = 0.8A. 

VIII Asymptotic behavior 

As a last issue, we will analyze the solutions of our model when R is very large in comparison with 
the de Sitter curvature Rds- This means that Eq. (|V.1[) can be approximated by 

^ cx)) ~ 7i?" , (VIII.l) 
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Figure 10: Plot of R/K for i?o = A. 

which is proved by the fact that a > 2 and, by setting 7 = -^ds^^, one has ^R" 3> R- In order to 
check for solutions, we use Eq. (|IL22[) and verify Eq. (jll.20[) . A class of asymptotic solutions of 
the model of Eq. (jVIILip at the limit t -J> 0+ is 

Hit)^^, (VIII.2) 

where Hq is a large positive constant and /3 a positive parameter so that f3 — 1 or f3 > 1. It follows 
from Eq. ((IlTqI) . 

i?^12||. (VIII.3) 



Eq. (flL22|) gets 



Poff ^ ^ . (VIII.4) 



Here, 5 is a positive constant and Eq. (jll.201) . in the limit t 0+, is perfectly consistent. This result 
shows that in the limit i? — 00 the model exhibits a past singularity, which could be identified 
with the Big Bang one. It is important to stress that this kind of solution is disconnected from 
the de Sitter inflationary solution, where the term R is of the same order of 7i?" and is therefore 
not negligible as in Eq. (jVIILip . We may assume that, just after the Big Bang, a Planck epoch 
takes over where physics is not described by GR and where quantum gravity effects are dominant. 
When the universe exits from the Planck epoch, its curvature is bound to be the characteristic 
curvature of inflation and the de Sitter solution takes over. 



IX On the stability of de Sitter space and a realistic model 
without singularities 

We here investigate in more detail the stability of the de Sitter solution (or its absence) and 
construct another model which does not generate any singularity. The de Sitter condition (III. lip 
can be rewritten as 

Let R = Rds be a solution of pX.ip . Then F{R) has the form 

m) 



R^ 



= /o + /i(i?)(i?-i?ds)" . (IX.2) 
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Here, /o is a constant, which should be positive if we require F{R) > 0. We now assume that n is 
an integer bigger or equal than 3; n > 3. Assume the function f{R) does not vanish at i? = i?dSj 
f{Rds) 7^ 0. Since n > 3, one finds 



i?dS 



F'jRds) 
F"iRds) 



0, 



(IX.3) 



which tells us that in pi.lSp vanishes. Therefore, a more detailed investigation is necessary in 
order to check stability. Using the expression of to^. 



R 



4F{R) 



F'iR) {F'{R)f F"{R) 



(IX.4) 



one can investigate the sign of in the region R ^ _Rds ■ Note that the Eq. (|II.15P is not used 
since this expression is only valid at the point R = Rds- Hence, we get 



3n(n-l)i?gs/i(j?ds) 



(IX.5) 



Eq. (|IX.5p indicates that, when n is an even integer, the de Sitter solution is stable provided 
fi{Rds) < but it is unstable if /i(i?ds) > 0. On the other hand, when n is an odd integer, the de 
Sitter solution is always unstable. Note, however, that when /i(-Rds) < (/i(i?ds) > 0), we find 
ml>0 {ml < 0) if i? > i?ds, but ml < {ml > O) if i? < i?ds- Therefore when /i(-Rds) <0, R 
becomes small but when /i(i?ds) > 0, R becomes large. The stability condition can thus be used 
to get realistic (unstable) de Sitter inflation for a specific F{R) gravity. 
The following model, instead of (|V.ip . is considered (compare with [23^), 



F{R) ^1 2A 
i?2 R r' 



aR- 



R^^ 
n 



R — Ri 
Ri 



1 



(IX.6) 



Here n is assumed to be an odd integer n > 3. We also assume that Ri, a, and e are positive 
constants and that e satisfies the condition < e < 1/n. We choose a to be small enough. When 
< i? <C we find i? ~ i? and therefore the model pV.ip is reproduced. When R ^ Ri, 
F{R)/R^ behaves as in p}C2|) . with 



Ri 



2nA 



Ri 



fi{R.) 



n 

Ri 



2nA 



R, 



Ri 
nR(] 



Ri 



Since < A ^ 1 and it is assumed that a <^ 



we find 



7% Ti 
Hi Hi 



(IX.7) 



(IX.; 



and, therefore, there exists a de Sitter solution R = Ri and the curvature always becomes smaller, 
slowly decreasing from the de Sitter point. Therefore, no future singularity is generated. When 
i? — >■ oo, F{R) behaves as 

F{R) - ai?2-2ne^e(«-l) ^j-j^ 

Since 1 < 2 — 2ne < 2, the singularity cannot emerge. Using the same numerical techniques as 
in the above sections, one can numerically fit this non-singular model with actual observable data 
coming from the dark energy epoch. 



X Discussion 



In summary, we have investigated in this paper some models corresponding to the quite simple 
exponential theory of modified F{R) gravity which are able to explain the early- and late-time 
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universe accelerations in a unified way. The viability conditions of the models have been carefully 
investigated and it has been demonstrated that the theory quite naturally complies with the local 
tests as well as with the observational bounds. Moreover, the inflationary era has been proven to 
be unstable and graceful exit from inflation has been established. A numerical investigation of the 
dark energy epoch shows that the theory is basically non-distinguishable from the latest observa- 
tional predictions of the standard ACDM model in this range. Special attention has been paid in 
the paper to the occurrence of finite-time future singularities in the theory under consideration. It 
has been shown that it is indeed protected against the appearance of such singularities. Moreover, 
its evolution turns out to be asymptotically de Sitter (it has a late-time de Sitter universe as an 
attractor of the system). Hence, the future of our universe, according to such modified gravity, is 
eternal acceleration. We have also demonstrated that slight modifications of the theory may lead 
to other non-singular exponential gravities with similar predictions, what points towards a sort of 
stable class of well-behaved theories. 

Very nice properties of exponential gravity are its extreme analytic simplicity, as well as the 
noted singularity avoidance. In this respect, the theory considered seems to be a very natural 
candidate for the study of cosmological perturbations and structure formation, which are among 
the most basic issues of evolutional cosmology. However, the theory remains in the class of higher- 
derivative gravities, which is not yet well understood, even concerning its canonical formulation 
|25j . In this respect, the covariant perturbation theory developed in [26j could presumably be 
applied for such investigation. This will be pursued elsewhere. 
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A The Einstein frame 



F{R) gravity may be rewritten in scalar-tensor or Einstein frame form. In this case, one can 
present the Jordan frame action of modified gravity of Eq. (jll.ip by introducing a scalar field 
which couples to the curvature. Of course, this is not exactly a physically-equivalent formulation, 
as explained in Ref. |27| . However, Einstein frame formulation may be used for getting some of 
the intermediate results in a simpler form (especially, when matter is not accounted for). 
Let us introduce the field A into Eq. (|II.1[) : 



^ / V^)[F'iA){R-A)+F{A)]d^x. 



(I.l) 



Here " Ji^" means "Jordan frame". By making the variation of the action with respect to A, we 
have A = R. The scalar field a is defined as 



a = -\n[F'{A)]. 

Consider now the following conformal transformation of the metric, 



(1.2) 



(1.3) 



for which Eq. p.ll) is invariant. By using Eq. (|I.3[) . we get the Einstein frame (EF) action of the 
scalar field a: 



I 



EF 



1 

2^ 



d X , 



(1.4) 
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where 

V(<7) = c^hic-") - c^" Flhic-")] . (1.5) 



/i(e"'^) is the solution of Eq. (O)) : 

h{e"')^A. (1.6) 

In order to pass to the scalar-tensor theory, we need the explicit form of the potential V{a). In 
principle, the result of Eq. (jI.6P will be in the form of a complicated transcendental function. 
However, in exponential gravity the calculation simplifies a lot. 

Thus, for the one-step model in Eq. pV.ip with A = i?o, Eq. ((L6| leads to 

/i(e-") = ln(l-e-")"'/^° , (1.7) 
and Eq. (11.51) is a simple transcendental equation which yields the nice result 

V{a) ^ - — e'"(l - e'^) In (l - e"'^) + 2Ae''(2i?o-i)/i?o _ (j g) 

Rq 
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